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A Helmholtz resonator can be used as a transmission part to connect the thermoacoustic engine (TE) with 
its load. However, the resonator can significantly influence the performance of the TE. In order to inves¬ 
tigate the impact of a Helmholtz resonator on the onset process of a TE, infrared (IR) imaging is firstly 
used as a visualization method to characterize the onset mechanism. The influence of dimensions of 
the Helmholtz resonator on the onset process are analyzed experimentally and theoretically. Results 
show that the Helmholtz resonator reduces the pressure amplitude at the onset moment and increases 
the onset temperature of the TE, both of which depend on the acoustic power absorbed from the TE. 
Onset without a sudden increase of pressure amplitude is observed with the Helmholtz resonator at res¬ 
onance length. This paper shows that IR imaging is an effective way to characterize the temperature dis¬ 
tribution in a TE study. 

© 2011 Elsevier Ltd. All rights reserved. 


1. Introduction 

A thermoacoustic engine (TE) with no moving parts is a new 
type of long-life, non-polluting engine that can utilize thermal en¬ 
ergy, even low grade thermal energy [1]. It can be used to drive 
coolers and generators in the fields of electronic device cooling, 
natural gas liquefaction, pulse tube refrigeration, and power gener¬ 
ation. Good performance of these cooling and power generation 
systems requires a well-designed transmission line to connect 
the thermoacoustic engine and the load or good phase shifter to 
improve the performance of the pulse tube cryocoolers (PTC). 
The transmission line or the phase shifter are commonly composed 
of a long slim tube, which has been investigated by researchers in 
the fields of PTCs and thermoacoustic systems. In 2005, an “acous¬ 
tical pump” [2] was proposed by Luo et al. to increase the pressure 
ratio from 1.25 to 1.47. A cascade acoustic amplifier [3] was also 
proposed by the present authors to further increase the output 
pressure ratio of acoustic amplifiers and significantly reduce power 
consumption. A Helmholtz resonator [4] was then studied by Sun 
to focus acoustic power in a TE. In addition, as good phase shifters, 
inertance tubes are widely used in high frequency PTCs [5,6], the¬ 
oretical and experimental research has been conducted to predict 
their characteristics [6-9]. 
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As tools to control steady, simple harmonic sound field that 
have tonal, narrow band spectrums [10,11], Helmholtz resonators 
have been exploited to enhance or attenuate sound fields since an¬ 
cient Greece [11,12]. They are widely used as mufflers in pipes and 
engines [13,14] to reduce noise and vibration. And the Helmholtz 
resonator in our previous study [4] was used as a transmission line 
to enhance acoustic wave, which amplifies the pressure amplitude 
from 0.20 MPa (p in ) to 0.32 MPa (p ou t), Pout/Pin = 1-60. Then another 
type of Helmholtz resonator composed of a long open tube, a res¬ 
ervoir and a long close-end tube (see Fig. 1) was studied by us, 
which has larger amplifying ability [15] (p in = 0.14 MPa, 
Pout = 0.45 MPa, Pout/Pin = 3.21). In these two Helmholtz resonators, 
the inertance of the long tube is mainly responsible for the pres¬ 
sure amplitude enhancement. 

However, the impact of a Helmholtz resonator on the onset pro¬ 
cess of a TE has not been well understood. Onset indicates the tran¬ 
sition of TEs from the stationary to periodic oscillating state when 
the heating temperature is larger than the critical value (the so- 
called onset temperature). As onset temperature determines the 
grade of the heat that a TE can utilize, the onset mechanism has 
been attracting more and more academic attention. The tempera¬ 
ture distribution of the regenerator by thermocouples is commonly 
characterized to analyze the onset mechanism in experimental re¬ 
searches. The distribution by this traditional way is not complete 
because the numbers of the thermocouples is always limited. 

In order to get the full-scale temperature distribution of the 
regenerator, infrared (IR) imaging is firstly used as a visualization 
method to characterize the onset mechanism in TE study. An IR 
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knowing any two of the variables p out , C/ out , Pin, and C/ in , the other 
two can be calculated. The transmission characteristics of the sys¬ 
tem can be described using the transmission matrix: 
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Fig. 1. Schematic diagram of Helmholtz resonator. 
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Fig. 2. The values of p 0 ut/Pm and input acoustic power dependence on L 3 . 

camera is used to measure the temperature distribution of the 
regenerator when the TE connects to the Helmholtz resonators 
with long tubes shown in Fig. 1. The transmission characteristics 
of the Helmholtz resonator are calculated to further explain the 
experimental phenomena. 


2. Theoretical analysis 


The method to simulate the transmission characteristics of the 
Helmholtz resonator with long tubes in thermoacoustic systems is 
detailed in our previous paper [16]. For a system coupling some 
ducts, the relations between the output pressure (p out ), volume 
flow rate (U out ) and inlet pressure (p in ), volume flow rate (U in ) 
can be expressed as: 
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where Ai, A 2 , A 3 , ..., or A n is the transmission matrix of each individ¬ 
ual duct. Once the transmission matrix is obtained, assume that by 
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where Z n is the series impedance per unit length, Y n is the shunt 
admittance per unit length, L n is the length of each tube, co is the 
angular frequency, p m is the mean density of the working gas, A is 
the cross-sectional area of each tube, a is sound speed, y is ratio 
of isobaric to isochoric specific heats, f v and f k are viscous function 
and thermal function, s s is boundary-layer approximation, and 
i = V^ T. 

For the Helmholtz resonator shown in Fig. 1, the transmission 
characteristics of the Helmholtz resonator can be obtained accord¬ 
ing to Eqs. (l)-(3). It should be noted that these equations are only 
applicable for laminar flow. However, the flow mentioned in this 
paper is turbulent. For turbulent flow, a turbulence correction 
should be made. So far, corrections for turbulence in an oscillating 
flow is not a mature science; there are only vague ideas for the cor¬ 
rections to turbulence. We proposed to make further corrections 
for turbulence when R/S v (R is the radius of the duct and S v is the 
viscous penetration depth of the gas) is large by increasing the vis¬ 
cous and thermal penetration depth, introducing a factor n, which 
is detailed in our paper [16]. 

In the simulations, the working gas of nitrogen is used with 
pressure of 2.2 MPa and frequency of 23 Hz. The diameters of tube 
Li and L 3 are 29 mm and 8 mm. The length of L\ is 1.08 m. The vol¬ 
ume of reservoir l 2 is 1 L, which can be considered as a short tube 
with an inner diameter of 99 mm and a length of 130 mm. The 
dimensions of L\ and L 2 are fixed and the length of L 3 changed. H ou t 
is 0 as the end of L 3 is blocked. Fig. 2 shows that the results of p out / 
Pin and acoustic power are dependent on the length of L 3 when the 
inlet pressure amplitude is 0.025 MPa. The value of p ou t/Pin in¬ 
creases with L 3 to the maximum, then decreases gradually. The 
length of L 3 corresponding to the peak value of p ou t/Pin is called res¬ 
onance length, which is 3.5 m. The calculation results indicate that 
the resonator absorbs the largest acoustic power at the resonance 
length. The calculation results agree well with the experiments 
(see the points in Fig. 2) under the same condition. The resonance 
length by the calculations is a little shorter than that by experi¬ 
ments, so the calculation results are always higher than the 
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Fig. 3. Schematic diagram of the traveling-wave TE with an infrared imaging camera. Thermocouples are located at Positions H and C, pressure sensor is located at 
Compliance (Position P). 










































































510 


B. Wang et al./Applied Acoustics 73 (2012) 508-513 



A :178s before onset B: 118s before onset 


C: 58s before onset 


D: onset moment 



E 60s after onset 


F: 118s after onset 


G: 182s after onset 


H: 240s after onset 


Fig. 4. IR images of the regenerator in the onset process. 


experimental data before the peak value, but lower than the exper¬ 
imental data after the peak value. 

3. Experimental results 

A traveling-wave TE is the TE in which the sound wave is trav¬ 
eling. It has much higher potential efficiency than the standing- 
wave TE (the other kind of TE). Fig. 3 shows the schematic diagram 
of the traveling-wave TE with an infrared imaging camera. The 
traveling-wave TE is detailed in the literature [17], which consists 
of a loop and a resonator tube. The loop is the main part converting 
thermal energy to the acoustic power, consisting of heater, cooler, 
regenerator, and phase adjusters such as the feedback tube and 
compliance. The pressure data acquisition system is used in the 
experiments, which consists of pressure sensors (linear silicon pie¬ 
zoelectric type, KPY 46R), a data acquisition clip (NIPC-1200,12 bit 
precision and a sampling frequency of 10 kHz), a computer, and a 
self-developed program for LabView 7.1 by National Instruments, 
Inc. The comprehensive error of the average pressure measure¬ 
ment is ±2.22E-02 MPa and that of the pressure amplitude mea¬ 
surement is ±6.43E-04 MPa. 

Heating power is adjusted by changing the charging voltage to 
the heaters, and is displayed by a digital dynamometer. With a 
heating power of 1000W, the uncertainty of the heating power 
measurement is ±11.5 W. 

An infrared imaging camera made by FLIR SYSTEMS is used to 
measure the temperature distribution of the regenerator. Visual 
images of the radiated heat of the objects can be obtained using 
the IR thermal imaging camera in the absence of light, thus the vi¬ 
sual temperature distribution can be measured. IR imaging method 
may be weak in measurement accuracy when parameters are not 
set appropriately, thus temperature by thermocouples is measured 
to compare with that by IR imaging. A K-type thermocouple is used 
to measure the heater temperature, and two thermocouples are 
also used to measure the temperature near the hot and cold end 
of the regenerator (labeled as Position H and Position C, shown 
in Fig. 3). 

The Helmholtz resonators with different L 3 were connected to 
the TE (with nitrogen as a working gas, where the pressure is 
2.2 MPa) and the temperature distribution of the regenerator was 
measured visually by an IR camera. The Helmholtz with 
L 3 = 5.7 m is chosen to illustrate the onset process. Fig. 4 shows 
the temperature distribution evolution of the regenerator in the 
onset process. K-sheathed thermocouples placed at Position C 


and Position H appear at the right side of each image. Using soft¬ 
ware ThermaCAM QuickReport 1.1, the temperature range and col¬ 
or of the images can be changed, and the temperature at any 
position of an image can be obtained. The display temperature of 
all the images in Fig. 4 is adjusted in the range of 20-200 °C, where 
the black region represents the lowest temperature, and the white 
shows the highest. Fig. 4D is the image at onset moment, while the 
images 3 min before (Fig. 4A-C) and 4 min after (Fig. 4E-H) the on¬ 
set moment are also shown. All the images display obvious color 
separation in the axial direction of the regenerator, indicating the 
existence of a temperature gradient along the regenerator, while 
the same color in the horizontal direction indicates the same tem¬ 
perature. With a constant heating power of 1000 W, the tempera¬ 
ture of the regenerator gradually increases. The TE onsets when the 
heater temperature reaches 168.1 °C, then the temperature of 




Time (s) 

Fig. 5. Temperature and pressure amplitude evolution in the onset process: (a) 
temperature comparison between thermocouples (TC) and the infrared camera (IR); 
(b) temperature evolution of the heater and regenerator; (c) pressure amplitude 
evolution at compliance. 








































B. Wang et al./Applied Acoustics 73 (2012) 508-513 


511 


Linear temperature distribution 



Fig. 6. The axial temperature distribution in the regenerator in onset process. 
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Fig. 7. The waveforms under different conditions. 


the regenerator continues to rise because of the non-stop heating 
power. 

Fig. 5 shows the dependence of simultaneous acquisition of 
pressure and temperature on time, where 0 s indicates the moment 
when the heater starts heating the TE. The time of A,B,..., H in 
Fig. 5 is consistent with that in Fig. 4. The temperature results 
are given in Fig. 5a (TC represents the thermocouples). The temper¬ 
ature data in the IR images at Position C and Position H are also col¬ 
lected and compared with the data by the thermocouples. At 
Position C, the results by the two methods coincide well. At posi¬ 
tion FI, they are close with a small deviation, especially at a higher 
temperature, as the temperature uncertainty of the thermocouples 
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Fig. 8. Onset waveform when L 3 = 3.7 m. 
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Fig. 9. The variation of pressure amplitude when L 3 = 3.7 m. 


increases with temperature. The results shows that measurement 
accuracy by IR imaging is acceptable. Fig. 5b and c gives the tem¬ 
perature and pressure evolution during the onset process. Before 
onset (A, B, and C) the temperature of the heater and regenerator 
gradually increases. Meanwhile the temperature difference be¬ 
tween the hot and cold end of the regenerator rises, reaching the 
critical value at the time of onset (D), where onset begins. 

The temperature distributions of the regenerator along the axial 
direction are collected (see the white line Fig. 4A) and displayed in 
Fig. 6 with the non-dimensional length L D as the horizontal axis 
and non-dimensional temperature T D as the longitudinal axis. 
The definition is L D = L/AL, where L is the distance away from the 
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cold end of the regenerator and A L is the whole length of the 
regenerator; T D = (T- T 0 )/(T H - To), where T, I H , and T 0 represent 
the regenerator temperature at the position of L , hot end and cold 
end. The temperature distribution of the regenerator is nonlinear 
before onset (Fig. 6A-C). This nonlinear distribution is mainly 
influenced by the large thermal resistance in the axial direction be¬ 
tween the stainless screens in the regenerator and large heat leak¬ 
age in the radial direction (in the radial direction thermal 
conduction of the stainless screen is good enough, in addition, 
the outside wall of the regenerator is not insulated). When the 
temperature difference between the hot and cold end of the regen¬ 
erator reaches the critical value, the working gas starts to oscillate. 
In a traveling-wave TE the penetration depth of the gas is much lar¬ 
ger than the characteristic dimension of the stainless screen, the 
sufficient heat transfer between them leads to the change of the 
temperature distribution after onset (see Fig. 6E—H), the tempera¬ 
ture distribution line goes closer and closer to the linear line. 

The waveforms at different length of L 3 are compared in Fig. 7, 
meanwhile the waveform when TE onset occurs without the Helm- 
holtz resonator is also shown. The display range of the pressure ra¬ 
tio is adjusted the same to represent the obvious difference under 
different conditions. The Flelmholtz resonator decreases the pres¬ 
sure amplitude of the TE when it onsets. The pressure amplitude 
jumps to the largest value without the Helmholtz resonator. As 
above mentioned analysis in Fig. 2, when the Helmholtz resonator 
is connected to a TE, the value of the pressure amplitude at the on¬ 
set moment depends on the acoustic power consumption by the 
resonator. The pressure amplitude has the largest value when 


L 3 = 0, as it needs the smallest acoustic power. After the TE onsets, 
the pressure amplitude decreases to a constant immediately after 
it jumps to the peak value because of the existence of the Helm¬ 
holtz resonator. The onset waveforms are nearly the same when 
L 3 = 2 m and L 3 = 5.7 m, as the acoustic power consumptions of 
these are very close; both are larger than that when L 3 = 0. The 
pressure amplitude comes to the smallest value (near 0 MPa) when 
L 3 = 3.7 m, as the TE onset near the resonance length of L 3 , the res¬ 
onator absorbs the largest acoustic power. The detail is repre¬ 
sented in Fig. 8. 

When L 3 is close to the resonance length, Fig. 8 indicates that 
the pressure amplitude increases smoothly without a jump. No 
obvious onset waveform can be distinguished in a longer duration 
of 300 s. This duration includes onset process which lasts about 
25 s. After onset, the pressure amplitude increases with the heating 
temperature, thus no steady state is observed. In Fig. 9, the smooth 
change of pressure amplitude also appears in the whole process, 
including both onset and damping, which differs from the sharp in¬ 
crease under the other conditions when L 3 is far from the reso¬ 
nance length (for example, an obvious jump occurs in Fig. 5c 
when the TE onsets). In the onset mechanism study, this phenom¬ 
enon that a TE onsets without a pressure amplitude jump has not 
been observed before in the literature. 

The onset temperature T H (the heater temperature at the onset 
moment, measured by thermocouples) changes with the length of 
L 3 , see Figs. 10 and 11. The Helmholtz resonator increases the onset 
temperature of the TE. Without the resonator, T H is the lowest. 
With the resonators, T H changes with the length of L 3 and its value 
is related to the acoustic power consumption. The onset tempera¬ 
ture increases in turn when L 3 = 0 m, L 3 = 5.7 m, L 3 = 2 m, and 
L 3 = 3.7 m. The reasons are similar with the above-mentioned dis¬ 
cussion of acoustic power consumption. T H is the largest when 
L 3 = 3.7 m, as the resonator absorbs the largest acoustic power 
from the TE, which delays the onset process of the TE. 

The IR images of the regenerator at the onset moment also show 
the obvious differences when the length of L 3 changes, as Fig. 11. 
The resonator near the resonance length has the largest effect on 
the performance of the TE. 

The Helmholtz resonator with long tubes has two functions: the 
first is to reduce the pressure amplitude of the TE, the second is to 
amplify the pressure amplitude within the resonator along its 
length. The first has negative and the second has positive effect 
on the pressure amplification effect, thus the pressure amplitude 
at the end of the resonator is a compromise of these two functions. 




Lj=0m, 7^=161.0*0 Ly= 2m, r H =182.7’C 



Z. 3 =3.7m, 7^=301.8*0 Without resonator, 7^=157.0*0 


Z. 3 =5.7m, 7^=168.1*0 


Fig. 11. IR images of the onset moment with different L 3 . 
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When vibration caused by thermoacoustic effect is not wanted, 
resonators should work at resonance length; In contrast, when 
amplification of pressure amplitude is preferred, resonator should 
work away from the resonance length. 

4. Conclusion 

Using an infrared imaging camera as a means of characterizing 
temperature distributions of thermoacoustic regenerators, the vi¬ 
sual IR images have been obtained for the first time in a TE. The 
new method can reveal the full-scale temperature distribution of 
the regenerator that the traditional method by thermocouples can¬ 
not. As the verification basis for the simulation, especially for a 
Computational Fluid Dynamics (CFD) study, this study supplies 
sufficient experimental data showing that the infrared imaging 
may be used as a powerful method in the TE study. 

The impacts of the Helmholtz resonator on the onset process of 
a TE are studied theoretically and experimentally. The resonator 
reduces the pressure amplitude at the onset moment, which de¬ 
pends on the acoustic power absorbed from the TE. The amplitude 
jumps to the largest value without the Helmholtz resonator. Mean¬ 
while, the Helmholtz resonator increases the onset temperature of 
the TE, which also depends on the acoustic power absorption. On¬ 
set without a sudden increase of pressure amplitude is firstly ob¬ 
served in the TE study when the Helmholtz resonator is at its 
resonance length. 
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